Synchronization in networks of slightly nonidentical elements 
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We study synchronization processes in networks of slightly non identical chaotic systems, for which 
a complete invariant synchronization manifold does not rigorously exist. We show and quantify how 
a slightly dispersed distribution in parameters can be properly modelled by a noise term affecting the 
stability of the synchronous invariant solution emerging for identical systems when the parameter 
is set at the mean value of the original distribution. 



PACS numbers: 
Keywords: 



.75.-k, 



'.75. He, 05.45.Xt 



Complex networks are the prominent candidates to de- 
scribe sophisticated collaborative dynamics in many ar- 
eas Recently, the dynamics of complex networks has 
been extensively investigated with regard to collective 
(synchronized) behaviors 0], with special emphasis on 
the interplay between complexity in the overall topology 
and local dynamical properties of the coupled units. The 
usual case considered so far is that of networks of iden- 
tical dynamical systems coupled by means of a complex 
wiring of connections. In this framework, several studies 
have shown how to enhance synchronization properties, 
by properly weighting the strengths of the connection 
wiring [J-Q. 

In this paper, we extend the study of synchronization 
phenomena in complex networks to the case of slightly 
non identical coupled dynamical systems, i.e. networks 
whose nodes are represented by dynamical systems each 
one of them evolving with the same functional form, but 
with a different, node dependent, value of the control pa- 
rameters. This is motivated by the fact that such a rep- 
resentation seems a more adequate description of many 
relevant phenomena occurring in natural systems, where 
the hypothesis that the evolution in different nodes be 
identical is very often a too restrictive assumption. 

Let us then consider a network of N coupled dynami- 
cal systems with slightly mismatched control parameters 
described by the equations 



N 

F(x ijgi )-a^G i . 



1,...,N, (1) 



where x, are the state vectors in each network node, F 
defines the vector field of the considered systems, gi are 
the control parameter vectors, H[-] is an output function, 
and a is the coupling strength. G is the Laplacian matrix 
of the network. As so, it is a symmetric zero row sum ma- 
trix, it has a real spectrum of eigenvalues Ai < • • • < A at, 
Gij (* 7^ i) i s equal to 1 whenever node i is connected 
with node j and otherwise, and Gu = — J2j^a 



When the considered network consists of identical ele- 
ments (i.e., gi = g, Vi) the stability of the synchronous 
state [xj(t) = x s (t), Vi] is known to be determined by 
the diagonalized linear stability equation Q, yielding N 
blocks of the form 



Ci= [JF(x s ,g)- C rAiJH(x s )]Ci ) 



(2) 



where J is the Jacobian operator. The blocks ([2]) differ 
from each other only by the eigenvalues Ai < ■ • ■ < Ajv 
of the coupling matrix G. Replacing a\ by v in equa- 
tion ((2J, the behavior of the largest (conditional) Lya- 
punov exponent A vs v [also called master stability func- 
tion completely accounts for linear stability of the 
synchronized manifold. Indeed, the synchronized state 
associated with Ai = is stable when all the remain- 
ing blocks related with the other eigenvalues Ai (i — 
2, . . . , N) of coupling matrix G arc characterized by the 
negative Lyapunov exponents. So, to analyze the sta- 
bility of the synchronized state x s (t) in the network ([1]) 
only one parametric variational equation 



C=[JF(x 5 ,g)-zA/H(x s )R 



(3) 



should be considered to obtain the dependence of the 
master stability function A on the parameter v. Fur- 
thermore, the vector state x s (i) may be obtained as a 
solution of the uncoupled equation 



x s (i) =F(x fl (i),g). 



(4) 



It is worth noticing that the master stability function 
A(v) may be negative for a finite interval of ^-parameter 
values I st = i v i] v 2) @ or for an infinite one (z^ 2 = oo), 
depending on the specific choice of the functions F and 
H. The stability condition is satisfied if the whole set 
of eigenvalues A,; (i = 2, . . . , N) multiplied by the same 
a falls into the stability interval I s t, i.e., when condi- 
tions 17A2 > V\ and ctAat < take place simultaneously. 
The vector functions F(-) and H[-] are determining the 
boundaries V\ and V2 of the stability interval I s t , while 
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FIG. 1: (Color online) The fragments of master stability func- 
tion corresponding to those v ranges for which A(^) crosses 
the horizontal axis. In both plots a) and b), the dashed line 
refers to A D (v) calculated for the network of Rossler oscilla- 
tors with slightly different parameters uii at D = 3.5 (see text 
for details), whereas the solid line depicts A(j/) calculated for 
the network of identical Rossler oscillators with Wi = Hi = 1. 
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FIG. 2: The dependence of the length L st of the stability 
interval I® t on the noise intensity D for the network of Rossler 
oscillators with slightly detuned parameter u)i. The value LfJ 
is normalized on the length L s t of the stability interval I st for 
the network consisting of the identical Rossler oscillators 



the eigenvalue distribution is solely ruled by the topology 
of the imposed wiring of connections. 

Natural systems, however, are modelled by networks 
that generally consist of elements for which parameters 
might differ. Therefore, equation §3§ cannot be seen as a 
suitable representation of this case. As soon as the vector 
gi depends on i, an invariant synchronization manifold 
Xj(f) = x s (t), Vi no longer exists, and therefore the argu- 
ments of the master stability function do not rigorously 
apply. However, it has been numerically verified in Ref. 
[3| that, when the difference in the parameters is limited 
to a slight mismatch, the synchronization behavior keeps 
on holding in the synchronization region predicted by the 
master stability function of the system corresponding to 
the parameter vector g = (gA, where (.) stays for the 
ensemble average on the network nodes. 

In the following, we will give ground to such a numeri- 
cal evidence, and show that, unless a rigorous treatment 
of the complete synchronization state is prevented, an 
approximate treatment of these synchronization phenom- 
ena is possible under the assumption of a smallness in the 
parameter mismatch. Without lack of generality, we will 
develop our points with respect to a subclass of chaotic 
systems, namely the class of functions F describing self 
sustained chaotic oscillators. 

The core idea that justifies our approximation comes 
from the well known property Q of a pair of coupled 
identical chaotic oscillators, where there are two impor- 
tant values o\ < 02 of the coupling strength a determin- 
ing the transition to complete synchronization. Precisely, 
<j\ determines the blowout bifurcation when the 

largest tangential Lyapunov exponent crosses zero. The 
second one, 02 corresponds to the loss of the stability 
in the tangential direction of the unstable periodic or- 
bits with the lowest period embedded into the synchro- 
nized chaotic manifold. When the coupling parameter 
value is in the interval <J\ < a < 02 the bubbling phe- 
nomenon ll|, [12( may be observed. If one considers two 



coupled identical oscillators the synchronous regime is 
detected (after expiration of the transient) for coupling 
strength values er > a\ . Alternatively, if the control pa- 
rameters of the coupled oscillators differ slightly from 
each other, the synchronous behavior may be detected 
for the coupling strength values exceeding the threshold 
CT2- The same effect takes place if two identical oscilla- 
tors in the presence of noise are considered Q- In both 
cases the onset of synchronization is shifted towards the 
larger values of coupling strength a and determined by 
the 02-value. 

The idea here is that the property for synchroniza- 
tion of the network consisting of elements with slightly 
mismatched control parameters gi may be estimated by 
means of the consideration of the network of identical el- 
ements with the control parameter g = (gA in the pres- 
ence of noise. In practice, this assumption means that 
one can still evaluate the conditional Lyapunov expo- 
nents by means of equation © at g = (gi), but the evo- 
lution of the state x s (i) around which the conditional 
exponents are evaluated has to be taken as a solution 
of a stochastic differential equation, i.e. the evolution 
equation Q has to be replaced by 



x,(t)=F(x.(f)) +!>£(*), 



(5) 



where D£(t) is a noise term with zero mean value. 

The new stability interval I® = (i / [ ) ,i / 2 D ) for the se- 
lected intensity D of noise may be found in the same way 
as it has been described above by means of calculating 
the master stability function K D (v). We will show that, 
under such an assumption, the increase of the noise in- 
tensity D leads the boundaries Vy and of the stability 
interval I® to converge to asymptotic values v\ and z/J, 
respectively. These points v\ 2 are analogous to the cou- 
pling strength value 02 in the case of two coupled chaotic 
oscillators and determine the stability interval I* t for the 
considered network of slightly non identical elements. 



3 



To illustrate the proposed approach let us consider 
a network of coupled Rossler systems. The dynam- 
ics of such network is described by equation ([1} with 
Xj = {xi,Vi,Zi), = (ui), F(xj,gj) = P(xj,w j ) = 
[—UiVi - Zi,uiiXi + 0.165j/i, 0.2 -I- Zi{xi - 10)], H[x] = x: 



.1 i 

Vi 

Zi 



■ ^^iVi — Zi — (7 X^"=l GijXj, 
LO l X l + 0.165j/ i; 

0.2 + Zi(xi - 10). 



(6) 



The mismatch in the parameters here corresponds to 
a detuning in the natural frequencies Ui of the oscilla- 
tors, that are supposed to be distributed randomly with 
a mean value ui = (wj) = 1 and a small dispersion 
Aw ^ 0.1. 

For the case of identical Rossler oscillators (i.e. as- 
suming all oscillators to have the same natural frequency 
Q = 1) , it exists a finite range of values for the parameter 
v (that we will call the stability interval I s t) for which 
the master stability function is negative [fjj]. 

To take into account the small difference in the fre- 
quencies of the coupled oscillators, equations (|3|) and ([5]) 
are instead used to calculate the stability interval I® as 
discussed above. To calculate the master stability func- 
tion h D {v) characterizing the property for synchroniza- 
tion of the network with slightly detuned elements, we 
made use of a random process £ (t) distributed uniformly 
over the interval ( — 1.0; 1.0). To inte grat e equation §5§ 
the one-step method has been applied [13[ with the time 
step At = 10~ 6 . 

The fragments of the dependence A D (v) around the 
boundary points V\$, are shown in Fig. [TJ One can see 
that the influence of noise in ([5]) results in the shift of the 
boundary points v^ 2 an d m the consequent reduction of 
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FIG. 3: (Color online) {E) (see text for definition) vs a for 
the network of Rossler oscillators with identical parameters 
uji = ui — 1 (black solid line) and for the network of Rossler 
oscillators with slightly detuned parameters Aw; = 0.07 (red 
dashed line with black circles). The arrows labelled with "1" 
delimit the interval of coupling parameter a for which the 
error goes to zero in the case of identical oscillators, while the 
arrows labelled with "2" delimit the interval of a for which the 
error goes to zero when the oscillators have slightly detuned 
parameters. 
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FIG. 4: The length Lf t of the stability interval 7^ vs. the 
maximum deviation of the control parameter value Aw for 
the network of Rossler oscillators with slightly detuned pa- 
rameters. The value Lf t is normalized to that of the length 
(L a t) of the stability interval I 3t for the network consisting of 
the identical Rossler oscillators. 



the stability interval I®. Therefore, the range of the cou- 
pling strength value a corresponding to the synchronous 
dynamics of the network of elements with slightly differ- 
ent values of parameters is less in comparison with the 
analogous network consisting of the identical elements. 

The dependence of the normalized length L®/L st 
(where = v 2 — and L st = v 2 — vi, respectively) 
of the stability interval I® on the intensity of noise D is 
shown in Fig. [2] One can see that, under the increase of 
the noise intensity D, the length of the stability interval 
converges to the value L® t * which does not depend 
practically on D-value. At the same time, the bound- 
aries Vi and i>2 of the stability interval I® t converge to 
the points v\ and i^, respectively. Therefore, the ob- 
tained interval I* t is the region of parameter ^-values 
corresponding the stable synchronized behavior of the 
considered network consisting of elements with slightly 
different parameter values. It is important to note, that 
the stability interval I* t is found when the noise intensity 
D is increased step-by-step. At the same time, if the D- 
value used in ([5]) is too large (e.g., the noise intensity is 
comparable to the amplitude of oscillations) , the dynam- 
ics of oscillator may be destroyed completely by noise, 
and, as a result, the boundary points of the stability re- 
gion I* t will not be detected correctly. In other words, 
there is a range of the reasonable values of the noise in- 
tensity D corresponding to the behavior of the network 
with slightly detuned elements. 

In order to show that the approximate solution of the 
master stability formalism is valid already for networks 
of relatively small size, a direct numerical simulation of 
Eq. ([B} has been carried out for various different val- 
ues of the control parameters. This calculation allows 
to find the boundary of the stability of the synchronous 
regime directly and to compare them with the analogous 
ones obtained before (based on Equations ((3]), (0). We 
consider here a network of TV = 5 non identical Rossler 
systems with Aw^ = 0.07. The coupling matrix G is 
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FIG. 5: The noise intensity D vs. the control parameter 
deviation Aw. 
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-3/ 



•1.6; A 3 



(7) 



-2.0; 



having as eigenvalues Ai 
A 4 w -4.0; A 5 « -4.4. 

In the simulations, the appearance of a synchronous 
state can be monitored by looking at the vanishing of the 
time average (over a window T) synchronization error 



(E) 



T(N — 



1 r t+T 



(8) 



In the present case, we adopt as vector norm ||x|| = \x\ + 
\y\ + \z\. Fig.[3]rcports the synchronization error versus a 
for a given network topology. The Figure comparatively 
reports the case of identical Rossler oscillators all with 
frequency w = 1, and the case of slightly non identical 
oscillators with frequencies distributed around the same 
mean w = 1 and with Aw = 0.07. One can see that 
the interval of a for which the error goes to zero reduces 
in this case, in accordance with the arguments extracted 
from the Master Stability Function description. 

A relevant issue concerns the possibility of establishing 
a quantitative correspondence between the noise inten- 
sity D in equation ([5]) and the dispersion of the control 
parameter values in ([5]). To clarify this point, we consider 
the ratio L®/L st between the length of the stability in- 
terval for non identical elements in the presence of noise 
and the same length for the case of identical systems. 
For the network of elements with the slightly noniden 
tical control parameters ©, the stability interval 
may be defined as the coupling strength range where the 
synchronization error JHJ vanishes. 

The dependence of the normalized length Lf t /L st of 
the stability interval I^j. with the value of the maximal 
deviation Aw is shown in Fig. 2] By comparison with 
Fig. O it is apparent how the curve in Fig. 2] is in 
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excellent agreement with its analogue depicting the de- 
pendence of the normalized length L®/L st on the noise 
intensity D. Notice that, as the the nonidentity of the 
network elements (Aw) increases, the length Lfj. strives 
for its asymptotic value Lf t * , which now does not de- 
pend on the specific Aw- value. Moreover, this limit value 
Lf t * is in the good concordance with L^ t * . Finally, in 
Fig. [5] the relationship between the noise intensity D and 
the control parameter deviation Aw is reported, showing 
that the network of non-identical oscillators can be suit- 
ably modelled by a noise addition to the synchronization 
manifold characterizing the evolution of the correspond- 
ing network of identical units. 

In conclusion, we have estimated the property for syn- 
chronization of networks consisting of equal elements 
with slightly different control parameter values. This 
study may be considered as the extension of the already 
known method of analysis of the behavior of the networks 
of identical elements. 
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